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A new technique is described by which light-induced gauge potentials allow systems of ultra-cold 
neutral atoms to behave like charged particles in a magnetic field. Here, atoms move in a uniform 
laser field with a spatially varying Zeeman shift and experience an effective magnetic field. This 
technique is applicable for atoms with two or more internal ground states. Finally, an explicit model 
of the system using a single-mode 2D Gross-Pitaevskii equation yields the expected vortex lattice. 



I. INTRODUCTION 

Condensed matter systems are replete with many- 
body effects, where the interactions between the in- 
numerable particles determine the basic physics of the 
system. Of late, ultracold atoms have demonstrated 
a range of basal condensed matter systems and ef- 
fects: Bose-Einstein condensation (BEC) E][2], Tonks- 
Girardeau gases [3j H] , the superfiuid to Mott-insulator 
transition [5], Berezinskii-Kosterlitz-Thouless physics [6 
in bosons; and the crossover from a Bose conden- 
sate to a Bardeen-Cooper-Schrieffer paired superfiuid in 
fcrmions [7]. Here, I discuss a technique to realize more 
complicated states where the charge neutral atoms be- 
have as charged particles in a magnetic field. 

In a magnetic field a 2D electron gas (2DEG) , can dis- 
play a range of exotic phenomena, including the inte- 
ger quantum Hall effect (IQHE) and fractional quantum 
Hall effect (FQHE). The FQHE states are exotic quan- 
tum liquids, where the lowest energy charged excitations 
are fractionally charged quasiparticles. While remark- 
able, ideas of charge and spin fractionalization are now 
well established in modern descriptions of strongly inter- 
acting quantum systems. More recently, exotic - non- 
Abelian - states useful for topological quantum compu- 
tation have been predicted, but remain experimentally 
elusive 8 . Still, experimental evidence now strongly sup- 
ports the existence of fractionally charged excitations in 
these systems [9], but evidence for their statistics is less 
conclusive [TUlIPl] . 

Systems of ultracold atoms are uniquely positioned to 
realize topological phases arising in strongly interacting 
2D systems, Fermi and Bose alike. The more exciting 
FQHE states are inevitably very delicate and only exist 
in the most clean systems, if at all. As with a 2DEG, the 
states of the system can be labeled by the filling factor 
v = the ratio of the 2D particle density n to mag- 
netic flux $ = eB/h. When c<la 2DEG can display 
a range of exotic phenomena, including various FQHE 
states. The primary challenge is to engineer a Hamilto- 
nian for which neutral atoms behave as charged particles 
in a magnetic field. 

This paper describes a new procedure for creating 
light-induced gauge potentials [12j [131 EH ESI EE], as 
a way to create an effective magnetic field. This work 
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FIG. 1: Geometry and level diagram to realize light-induced 
gauge potentials, a. The figure depicts the two counter- 
propagating laser beams with momentum ki and k2, aligned 
along +x and —x respectively, and the physical magnetic field 
along y. b. Level diagram for two electronic ground states, 
c. Level diagram for three electronic ground states; a small 
quadradic Zeeman shift e of state |2) is not depicted. Addi- 
tional transitions, relevant for physical atoms, contribute to 
the state-independent light shift in this configuration. 



focuses on the adiabatic eigenstates of atoms in the pres- 
ence of two optical fields, and finds that the resulting 
Hamiltonian can describe charged particles in a mag- 
netic field. In contrast to earlier proposals, the requi- 
site spatial inhomogeneity here is provided by an exter- 
nal magnetic field gradient instead of inhomogenous op- 
tical fields. This paper first presents explicit results for 
a model system with two coupled states, and then gen- 
eralizes to the three-level case (relevant to the F = 1 
manifold of 87 Rb). The effective magnetic field result- 
ing from light-induced gauge potentials exists in a small 
spatial region. An important finding is that coupling be- 
tween more than two internal states increases the spatial 
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range over which large effective fields can be realized. 

In second-quantized form, the Hamiltonian for a parti- 
cle in a uniform time-independent magnetic field normal 
to a 2D plane is 
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where i/^( x ) is the field operator for the creation of a 
particle at x, and k xy — —id xy . For real magnetic fields 
the vector potential has gauge freedom. For example the 
Landau gauge choice, {A x — By,A y = 0}, gives a uni- 
form magnetic field along z. This proposal explicitly re- 
alizes Eq. [I] in a specific gauge (the Landau gauge for 
the geometry discussed below) : a Hamiltonian where the 
minimum of the energy-momentum dispersion relation 
E(k) becomes asymmetric [17] and is displaced from zero 
momentum as a function of spatial position. The dressed 
single-particle states are spin/momentum superpositions 
whose state decomposition depends on the local value of 
the effective vector potential A. In this way, the canoni- 
cal momentum associated with the Landau gauge is phys- 
ically observable by probing the internal-state decom- 
position of the adiabatic dressed states. This effective 
Landau-gauge vector potential was recently measured by 
Lin et al [15]. 

Our approach relies on a collection of bosons with two 
or more relevant electronic ground states interacting with 
two counter-propagating "Raman" lasers aligned along x 
that are detuned from each other by Ar, shown in Fig. 
[l] A small magnetic field B = (Bo + AB)y, introduces a 
linear Zeeman splitting between the levels; Ar — g^B^, 
so 5 = g^AB is the detuning from Raman rasonance. 
Here g is the atomic ^-factor, and [Lb is the Bohr mag- 
neton. I focus on the limit when both Raman beams are 
far detuned from the ground to excited state transition 
so there is negligible population in the excited state, and 
the Raman beams simply induce a coupling ft between 
ground states. As shown below, this set of coupling fields 
can lead to effective magnetic fields. 



II. TWO-LEVEL SYSTEM 

For simplicity, first consider a two level system with 
internal states |+) and |— ), where exact solutions, stud- 
ied in the context of a 3D BEC in Ref. [17], are read- 
ily available. (Physically, these two states might be 
two mp levels in the ground state manifold of an alkali 
atom; for example, the F = 1 manifold of 87 Rb at large 
enough field that the quadratic Zeeman effect resolves 
the three Zeeman sublevels.) Since counter-propagating 
Raman beams aligned along x couple states differing in 
k x by 2k r , the recoil momentum k r — 2t:/X and energy 
E r = h 2 k r 2 /2m will be taken as the units of momentum 
and energy. Here, A is the wavelength of the nearly de- 
generate Raman beams, m is the atomic mass, and the 



two- photon Raman coupling is fi. In the frame rotat- 
ing at Aft/h the Raman fields are detuned S = g^sAB 
from resonance and the atom-light coupling term in the 
rotating wave approximation (RWA) is 
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The notation ^(/c^, y) denotes the creation of a particle 
with wave vector k x along x at position y, with a = 
± and H.C. indicates the Hcrmitian conjugate. Also, 
observe that H' includes the Raman detuning terms. In 
the following analysis ft and 5 will be treated as spatially 
varying functions of y, but not x. 

Absent coupling, the Hamiltonian for particles in 2D 
is a sum H = H x + H y + V + H int . Respectively, these 
represent motion along x, motion along y, the external 
potential, and interparticle interactions. When expressed 
in terms of the real space field operators VvM, these 
terms are 
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The contact interaction for collisions between ultra-cold 
atoms in 3D is set by the 3D s-wave scattering length 
a s , here assumed to be state independent. Strong con- 
finement in one direction yields an effective 2D coupling 
constant g2D = V8nh 2 a s /mliio- lno is the harmonic 
oscillator length resulting from a strongly confining po- 
tential along z; a ID optical lattice, for example [TS]. Fi- 
nally, V(r) is an external trapping potential, also taken 
to be state-independent. 

This problem is exactly tractable when considering free 
motion along x, i.e., treating only H x and H'. The sec- 
ond quantized Hamiltonian for these two contributions 
can be compactly expressed in terms of the operators 
{0l(k x ,y),0l(k x ,y)} = {$ + (k x - l,y),4>l(k x + 1,2/)}. 
Using this Nambu spinor, H « H x + H' reduces to a 
integral over 2x2 blocks 
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labeled by k x and y. The dependence of the two-photon 
coupling ft and detuning 5 on y has been suppressed for 
notational clarity. The resulting Hamiltonian density for 
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motion along x at a fixed y is 

For each k x , H(k x ,y) can be simply diagonalized into 
spin-momentum superposition states by the unitary 
transformation U(k x ,y)H(k x ,y)U'(k x ,y)- The resulting 

eigenvalues E ± (k x ,y) = k 2 x + 1 ± \J (4k x - S) 2 + ft 2 / 2 
give the effective dispersion relations in the dressed basis 
$ a {k) — Uca'ity&a'ik)- Such states have been exten- 
sively studied in the context of velocity selective coherent 
population trapping, and for each k the eigenvectors of 
H(k x ) are said to form a family of states 20J. In terms 
of the associated real-space operators tp'^ (r) these diago- 
nalized terms of the initial Hamiltonian are 



H x + H' 
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In analogy with the terms "band" and "crystal- 
momentum" for particles in a lattice potential, the set 
of states giving rise to each dispersion curve will be 
called a "quasi-band" , and the quantum number k x the 
"quasi-momentum" . Here —id x is the real-space rep- 
resentation of the quasi-momentum k x . The symbol 
E± (—id x ,y) is a differential operator describing the dis- 
persion of the dressed eigenstates, just as the operator 
E x (—id x ,y) = (—ihd x — eBy) 2 /2m describes quadratic 
dispersion along x of a charged particle in the Landau 
gauge. 

To lowest order in 1/ft and second order in k x , 
E±(k x ,y) can be evocatively expanded: 
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Atoms in the dressed potential are significantly changed 
in three ways: (1) the energies of the dressed state atoms 
are shifted by e± = 1 ± ft/2 (the additional energy offset 
of order S 2 is relevant to trapping); (2) atoms acquire an 
effective mass m* /m = fl/(ft ± 4), and crucially (3) the 
center of the dispersion relation is shifted to eA x /hk r = 
J/(f2±4). A x can depend on y by virtue either of a spatial 
dependence on ft (as in Refs. [UlEI], which focused on 
the large ft limit), or via 5(y) — gfiB^B(y) as described 
below. In either case, the effective Hamiltonian is that 
of a charged particle in a magnetic field expressed in the 
Landau gauge. 

Figure [2] shows the dressed state dispersion relations 
from this model. Panels a and c show the undressed case 
(ft = 0) for detuning 5 — and 5E r respectively. Panels 
b and d depict the same detunings, for ft — l6E r , where 
the exact results (solid line) are displayed along with the 
approximate dispersion (red dashed line). Fig. [2Jd then 
shows the strongly dressed states for large ft, each of 
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FIG. 2: (Color online) Each panel pictures the dressed state 
dispersion relations for a two level dressed state atoms. The 
horizontal axis is quasi-momentum k x , and the vertical axis is 
the dressed state energy E±(k x ). The black lines are the exact 
eigenvalues of Eq. [2] and the red dashed lines are the analytic 
approximation, (a) bare potentials (undressed) with Raman 
beams on resonance; (b) dressed potentials (fl = 16E r , 5 = 
0E r ); (c) bare potentials (undressed) with Raman beams off 
resonance (S = 5E r ); and (d) dressed potentials with Raman 
beams off resonance (fl = 16E r , 8 = 5E r ). 



which is symmetric about k = 0, when detuned as in 
panel c, the dispersion is displaced from k — 0; when 
spatially dependent this displacement leads to a non- 
trivial gauge potential. In the limit of very small ft the 
dressed curve E_(k x ,y) forms a double-well "potential" 
as a function of k x ; in a related Raman-coupled system, 
Bose condensation in such double-well potentials have 
studied theoretically in Refs. [T71 |2"TI [2^1 12"3] . 

It is also possible to treat this problem in the Born- 
Oppenheimer (BO) approximation in which only H' is 
diagonalized [TJJ [T3J [TJ] [24"] . In this new eigen-basis, 
H x has off-diagional terms which are taken to be small, 
and ignored in the BO approximation. Such an assump- 
tion is valid only when ft~ l is small, in which case the 
BO approximation yields a dispersion exactly in the form 
of Eq. [4] where the crucial terms are mj /m = 1 and 
eA BO /hk R — ±5/^/5 2 + ft 2 . These relations converge to 
those in Eq. [4] for very small ft^ 1 . In contrast, the BO 
approach fails to yield the correct physics in the limit of 
small coupling, for example never predicting the double- 
well structure in E-(k x ). 



A. Effective fields, trapping, and optimization 

The analysis of the Raman coupling lead us to a 
dressed dispersion along x, and as is shown below, mo- 
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tion along y is largely unaffected. When the detuning 
is made to vary linearly along y, 6(y) = 5'y, an effec- 
tive single particle hamiltonian contains a 2D effective 
vector potential qA/hk r w {5'y /{A ± O), 0} - the vector 
potential for a magnetic field normal to the x-y plane ex- 
pressed in the Landau gauge. The effective magnetic field 
is q(3 z /hk r w (S'/(4± f2). Figure [3] shows the computed 
vector potential as a function of detuning 8. As expected, 
the linear approximation discussed above (dashed line) is 
only valid for small 5; as a consequence the effective field 
decreases from its peak value as S increases (top inset). 

In addition this technique modifies the trapping po- 
tential along y, i.e, it produces an (unwanted) scalar po- 
tential in addition to the vector potential. When the 
initial potential V(x,y) is harmonic with trapping fre- 
quencies lo x and uiy, the combined potential along y be- 
comes V±(y) = m(uj y + (uo* ± ) 2 )y 2 /2, where m(w±) 2 /2 = 
<S' 2 (4±ft)/4(4 + ft) 2 . 

This contribution to the overall trapping potential is 
not unlike the centripetal term which appears in a ro- 
tating frame of reference, where an effective magnetic 
field /3 ro t arrises as well. In the case of a frame rotat- 
ing with angular frequency f2 ro tj the centripetal term 
gives rise to a repulsive harmonic term with frequency 
w ro t — (<?A-ot/2m). In the present case the scalar trap- 
ping frequency can be rewritten in a similar form lu*± = 

(q(3/2m) x |4 ± fi| 3/2 /(4 + Q); the scalar potential may 
be attractive or repulsive, and it increases in relative im- 
portance with increasing fi. 

The effective field generated is inhomogeneous, how- 
ever, for many physical effects in a magnetic field, such as 
the Hall effect (quantum and classical) the filling fraction 
v, not the magnetic field, is the most relevant parame- 
ter. In the Thomas-Fermi limit the spatial density n of 
a BEC decreases quadratically from the center of a har- 
monic trap, at lowest order this can compensate for the 
decreasing effective field, leading to an extended region of 
constant v = hnjefiva. the systems center (bottom inset 
of Fig. 3j. This approach is very well suited for incom- 
pressible QHE states which will form a shell-structure at 
constant v; thus the effective homogeneity is enhanced 
for a harmonically trapped gas in the presence of a field. 
Still, this requires fine-tuning of the system-size for every 
/?; in the three level case, this fine-tuning restriction is 
lifted. 



fi = 16£ r 



B. Additional Coupling 



The preceding calculation omitted motion along y, in- 
teractions, and an external potential: H y , H lnt , and V. 
These three complicating terms can be treated easily, 
starting with the state-independent trapping potential 
V . The potential can be expressed in terms of dressed 
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FIG. 3: Effective vector potential eA x /hk r versus detuning 
5 for Q = 16E r . The solid line is the exact result and the 
dashed line is the lowest order expansion in 5. The insets 
depict predicted quantities relevant to experiment, computed 
for 87 Rb with a detuning gradient 5'(y) = 900 Hz//im, and 
A = 800 nm Raman lasers. The top inset shows the effective 
magnetic field B, indicating the degree of field inhomogeneity. 
The bottom inset shows the filling fraction v for a Rtf = 
37 (iim BEC chosen so v is nearly constant near the systems 
center. 



quasi-momentum operators via the relations 
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where V(Sk x ,y) is the potential Fourier transformed 
along x. The dressed states experience the same state 
independent potential as the initial states, however, 
the small off-diagonal terms of U(k\,y)W (k 2 ,y) to- 
gether with V{k\ — k 2 ,y) give transition matrix elements 
2V{k\ — k 2l y){ki — k 2 )/Vt between dressed states (through 
second order in f2 -1 ). In real space this gives a coupling 
oc il~ 1 d x V(x, y), sensible because the coupling term re- 
sults from deviations from a uniform potential. Since a 
typical trap is many tens of wavelengths in extent, and 
ultracold atoms generally have momenta at or below k r 
this coupling term is small, but not in general negligible. 

The term describing motion along y also leads to cou- 
pling terms between dressed states. The argument lead- 
ing to Eq. [5] for H y gives 



H y *-J d 2 r£^t (r) ^ (r); 



(6) 
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again having made the approximation U{k,y)dyW (k,y) 
w 0; this gives coupling (20) — 1 (d y S(y))d y at lowest order 
in This term results from a breakdown of a BO ap- 
proximation implicit in the diagonalization of H x + H' at 
fixed y leading to Eq. depends on y through 6(y)). 

This approximation is distinct from the BO approxima- 
tion alluded to earlier where the coupling Hamiltonian 
H' alone was diagonalized at fixed x and y. 

Finally, the arguments given above also show that the 
interaction H lnt leads to a dressed-state independent in- 
teraction with the same g with a state-changing coupling 
term proportional to c/2d/^ at order Q . 

Together these allow the construction of the expected 
"real space" Hamiltonian in the basis of localized spin- 
superposition states V>'(r). The BO violating coupling- 
terms discussed above can be treated perturbatively 
for particles in the lowest quasi-band leading to sta- 
ble cigenstates, however, transitions for particles start- 
ing in higher quasi-bands are energetically allowed and a 
Fermi's Golden Rule argument thus gives rise to "decay" 
from all but the lowest energy dressed state [2"5] . 

Using the standard argument of a single macroscopi- 
cally occupied state, the Hamiltonian reduces to the 2D 
Gross-Pit aevskii equation (GPE) 

{[E(-id x ,y)-%]+ V(r)+ 

(N - l) 32D |vI/(r) | 2 }*(!■) = /i*(r), 

using the x dispersion E(—id x ,y) which parametrically 
depends on y from Eq. [3j Thus one expects the usual 
formation of a vortex lattice at small effective fields when 
the single mode approximation is valid. 



C. Limitations 

Naturally, this technique is not without its limitations. 
Foremost among them is the range of possible qA x /hk r 
shown in Fig. [3] where = l6E r : while the linear expan- 
sion (dashed) is unbounded, the exact vector potential is 
bounded by ±1. The reason for this is clear; for example, 
the hybridized combination of |+, k — 1) and \ — ,k + 1) 
cannot give rise to dressed states with minima more pos- 
itive than k = +1 (where |1) is minimized absent dress- 
ing, Fig. |2j;); nor can the minima be more negative than 
k = -1. 

This limitation does not effect the maximum attainable 
field, only the spatial range over which this field exists. 
Specifically, a linear gradient in S(y) gives rise to the ef- 
fective field (3 z {y) which is subject to q(3 z (y)dy = 
2hk r . This simply states that the vector potential - 
bounded by ±hk r /q - is the integral of the magnetic field. 
Note however, that along x the region of large j3 z has no 
spatial bounds. 

A second limitation of this technique is the assump- 
tion of strong Raman coupling between the Zeeman split 
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FIG. 4: (Color online) Each panel denotes the dressed 
state dispersion relations E(k x ) for atoms dressed by counter- 
propagating Raman beams. The horizontal axis is quasi- 
momentum k x , and the vertical axis is energy in the RWA. 
(A) bare potentials (undressed) with Raman beams on res- 
onance; (B) dressed potentials (Q = 32E r , 5 = 0-E r ); (C) 
bare potentials (undressed) with Raman beams off resonance 
(5 — 5E r ); and (D) dressed potentials with Raman beams off 
resonance (fi = 32E r , 5 — WE r ). 

states. In the alkalis, when the detuning from atomic res- 
onance is large compared to the excited state fine struc- 
ture the two-photon Raman coupling for Amp = ±1 
transitions drops as f2 oc A~ 2 , not A -1 as for the AC 
Stark shift. As a result, the balance between off-resonant 
scattering and fl is bounded, and cannot be improved by 
large detuning. While this is a modest problem for rubid- 
ium (15 nm fine structure splitting), it is extremely im- 
portant for atoms with smaller fine structure splittings: 
potassium (sa 4 nm) and lithium (as 0.02 nm). This issue 
can be avoided for the two- level case, by using Am/ = 
transitions, e.g., between ground state hyperfine mani- 
folds in the alkalis. 



III. THREE-LEVEL SYSTEM 

The range of possible effective vector potentials can 
be extended by coupling more states, for example the 
mp states of a F > 1/2 manifold in the linear Zeeman 
regime. The calculation follows the two-level example 
above, except for the lack of compact closed- form solu- 
tions. Additional levels extend the range of the vector 
potential from ±fc r in the two level case to ±2Fk r for 
arbitrary F. 

For specificity, consider an optically-trapped system of 
87 Rb atoms in the F = 1 mainfold in a small magnetic 
field which splits the three mp levels by g/is|B| (Fig. [I]:). 
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FIG. 5: The left panel shows the optimized effective fields /3 
and the right panel the optimized filling fraction v. In each 
plot, fi = 16E r and the solid curve shows the case optimized 
for uniform field (e = —1.657E r ) and the dashed curve for 
uniform v (e = —0.064). 



The coupling fields can be produced by a pair of far- 
detuned counter-propagating lasers (aligned normal to 
the bias field B) detuned from each other lo\ — — 
<7/xb|B|//i — S. Laser polarizations, (cr+ + er_)/\/2 and 
7r, allow Raman transitions between the hyperfinc levels 
when the detuning A from the excited states is compara- 
ble or smaller than the 15 nm excited-state fine-structure 
splitting. 

As with the two level case, the ID Hamiltonian describ- 
ing motion parallel to the dressing lasers can be made 
block-diagonal. The 3x3 blocks H{k x ) describing the 
three internal states of the F — 1 manifold are 

/(k x -2) 2 + 5 ft/2 \ 

H{k x ) = ft/2 k 2 x + e ft/2 . (7) 

V ft/2 (k x + 2) 2 - S J 

In this expression, S is the detuning of the two photon 
dressing transition from resonance; e accounts for any 
quadratic Zceman shift; ft is the two-photon transition 
matrix element; and k x , in units of the recoil momen- 
tum k r , is the atomic momentum displaced by a state- 
dependent term k x = k — 2 for mj? = —1, k x — k for 
772 p = 0, and k x = k + 2 for m F = +1. When ft > 8\/2 
the three eigenvalues, denoted by E± and E are approx- 
imately: 

E ~( n V 1 (k V 2 ^ ±n 

± ~ \n±8V2j \ x ~ 8V2±n J + Vl 

(8) 

E « k 2 + 4. 

As with the two-level case, the states associated with 
eigenvalues E± experience an effective vector potential 
which can be made position-dependent with a spatially 
varying detuning 5 [26 . Again, a magnetic field gradient 
along y gives S cx y, and generates a uniform effective 
magnetic field normal to plane spanned by the dressing 
lasers and real magnetic field B. 



The resulting magnetic field is inhomogeneous and de- 
parts quadratically from its peak value. For experiments 
requiring constant filling fraction, v can be made ap- 
proximately uniform by proper selection of the system's 
Thomas- Fermi radius Rtf- Still, some experiments do 
require a homogenous effective field. In the lowest en- 
ergy quasi-band, the proper choice of e* = 4 - \/2ft/4 
(exact) suppresses the drop-off of the effective magnetic 
field, leaving terms of order <5 4 and higher (e results from 
quadratic Zeeman shifts and is controlled by the bias 
magnetic field B). For large enough ft, e* < corre- 
sponding to the physical sign of e in the 87 Rb F = 1 
manifold. For experiments benefiting from constant v, 
a similar analysis shows the filling fraction can be made 
constant to 0(5 4 ) when when e* = 3.2 — 0.204ft and when 
the usual Thomas-Fermi density profile goes to zero at 
<5 = 1.00ft — 8.5 (Fig. [5] shows an example of this opti- 
mization as well). 



A. Gross-Pitaevskii Equation 

The arguments leading to the GPE equation in the 
two level case remain valid here, and the coupling terms 
remain of the same order. A numerical solution to the 
GPE in the low- field regime is shown in Fig. [6j This cal- 
culation was performed for the lowest-energy of the three 
dressed states using the exact dispersion resulting from 
the numerical diagonalization of Eq. [7] The computa- 
tion uses 87 Rb parameters and the experimentally real- 
istic ft = 16E r . The inset to Fig. [6] depicts a case with 
trapping frequencies oj x /2tt = 10 Hz and Lo y /2ir = 40 Hz 
(the asymmetry of these terms is partially counteracted 
by the effective anti-trapping term along y resulting from 
the zero-offset of the dress state dispersion). The com- 
puted vortex lattice explicitly demonstrates that the ap- 
proach described above creates an effective field for neu- 
tral atoms in a non-rotating frame, even given realis- 
tic parameters. The main panel plots the inverse vor- 
tex spacing as a function of gradient directly obtained 
from the 2D GPE solution (symbols); overlapping these 
points is a solid line depicting the expected vortex spac- 
ing at the systems center (peak effective field) obtained 
by direct diagonalization of Eq. [7] The dashed line is 
the approximate expression from Eq. [8] The formation 
of the vortex lattice with the correct spacing clearly in- 
dicates that this technique does give rise to the expected 
effective magnetic field. 

A counterintuitive reminder of this simulation is 
that spatially stationary solutions to the dressed-state 
many body problem exist even when the lowest quasi- 
band wave-functions intrinsically involve large momen- 
tum components and spin-mixtures. To understand this 
situation we can consider a more pedestrian example: 
atoms in an optical lattice. In this case the systems sin- 
gle particle eigenstates - Bloch states - involve only one 
spin component but are composed of many momentum 
components each separated by 2k r . In the lowest band 
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FIG. 6: Inverse vortex spacing versus detuning gradient S'(y) 
in a three level system. The displayed symbols are obtained 
by solving the GPE for 3D BEC with 3 x 10 5 to 5 x 10 5 in the 
presence of an effective magnetic field with a detuning 5 — at 
the system's center, and a coupling f2 = 16E r . The simulation 
assumes a 10 pm Thomas- Fermi radius along z, and solves a 
2D GPE along the remaining two directions. The typical 
vortex spacing is obtained from the Fourier transform of the 
density distribution ^(r)] 2 . The uncertainties, reflecting the 
vortex-spacing distribution, are obtained from the half-width 
of the first peak in the same Fourier transform. The dashed 
line is the approximation, Eq. [8j and the solid line results 
from numerical diagonalization of Eq. [7| Inset: Calculated 
in-situ density distribution for a detuning gradient 5'{y) = 
0.023 x E r k r showing the expected vortex lattice structure in 
a non-symmetric and non-rotating system. 



of a sinusoidal lattice the q — Bloch state has no cen- 
ter of mass motion, and instead its many momentum 
components combine to produce the spatially periodic 
density modulation characteristic of Bloch states. In the 
present case of spin-momentum dressed states, the dif- 
fering momentum components neither result in center of 
mass motion, nor in density modulations as with an op- 
tical lattice. Instead, the momentum is associated with 
a spatially modulated spin texture aligned along x. In 
both cases, states away from local minima, with non-zero 
group velocity, do have non-zero mechanical momentum. 
The current case differ from the lattice analogy in one 



substantial way: here the analysis was performed in a 
frame rotating at the frequency difference between the 
Raman beams An/h. In the rotating frame the spin 
texture is static, however, in terms of the bare-states 
the time-dependent phase factors exjp(imF^Rt/h) imply 
that the local orientation of spin texture is rapidly vary- 
ing. 

IV. CONCLUSIONS 

Neutral atoms in the presence of suitable coupling laser 
fields experience effective magnetic fields and the explic- 
itly calculated coupling terms between dressed states are 
negligible only for atoms in the lowest energy dressed 
state. The same effective field effect exists for systems 
with three or more levels. This both extends the appli- 
cability of the technique, but in addition the additional 
level increase the spatial extent over which high effective 
fields can be realized. 

I discussed two optimizations: (1) where either the 
two- or three- level system can be fine-tuned to make the 
filling fraction constant through third order in displace- 
ment along y; and (2) where the quadratic Zeeman term 
in the three-level case allowed the field to be made uni- 
form through third order in y (the same type of reasoning 
could make v constant through order-5 in y) |27j . 

Finally, I showed that an explicit solution to the GPE 
equation in the presence of an effective gauge potential 
has the expected vortex lattice. Using this technique it 
is possible to generate effective magnetic fields sufficient 
to enter the FQHE regime, v < 1 where the GPE is 
invalid. For example, in the three level case optimized 
for uniform field (using Rubidium parameters, f2 = 16£" r 
and e = — 1.656£V as in Fig. ^ a gradient of 5'{y) » 
h x 16 kHz//im s» i.5E r /fjxa, requires a modest labora- 
tory magnetic field gradient of 2.3 T/m. This yields an 
effective field (3 = 4.14 mT, where the "effective charge" 
was taken to be e (in our technique the product qj3 is 
defined), and a magnetic length Ib — \Jhjqfi = 0.5 \xm. 
With a reasonable 2D atom density n as 1 /im~ 2 the fill- 
ing fraction is v = hn/q/3 — 1. Thus this approach allows 
experiments to reach the strongly correlated regime with 
realistic experimental parameters. 
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